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Abstract 

The variances of the quantum-mechanical noise in a two-input-port Michelson interfer- 
ometer within the framework of the Loudon-Ni model were solved exactly in two general 
cases: (i) one coherent state input and one squeezed state input, and (ii) two photon number 
states inputs. Low intensity limit, exponential decaying signal and the noise due to mixing 
were discussed briefly. 


1 Introduction 

In 1981 the effects of intensity fluctuations in the two light beams and radiation pressure on the 
mirrors in a Michelson interferometer were modeled in a unified way by Loudon [1]. In 1987, Ni 
extended the model to include the intrinsic uncertainties of the mirrors and obtained an exact 
solution of the variance of the quantum mechanical noise for a coherent light source with arbitrary 
intensity [2], These results were used recently by Ni in proposing an experimental scheme for 
controlling a macroscopic quantum state of a mirror by light shining [3], 

Quantum-mechanical noise of a Michelson interferometer is an important noise source in gravi- 
tational waves detection. Experiments had reached the shot noise limit already. Photon shot noise 
decreases as the intensity goes up. Roughly speaking, it is propotional to the inverse square root 
of intensity. But, it was argued [4] that the fluctuation of radiation pressure on the mirrors would 
increase as the intesity becomes higher so that a minimum would be reached, called the “standard 
quantum limit” [4]. On the other hand, various works [5] show that such a measurement (without 
loss) implies no limit while squeezed states were used. 

With the knowledge of squeezed states and that semi-classical model of interaction between 
a macroscopic object and photons, we got chance to probe this problem in detail. Within the 
framework of the Loudon-Ni model, the method of Ref. [2] was extended to obtain an exact 
solution for the variance of a two-input-port Michelson interferometer where a squeezed-state 
light source and a coherent light source (both with arbitrary intensity) were applied on each port 
respectively. Faithful matrix representations [6] were used in this calculation. The final result 
is more complicated. Nevertheless, it can be organized in a similar form as in Ref. [2]. Photon 
number state inputs can be treated in a similar way. 

Due to the complicated results we’ve got, the physical implication is yet under study. However, 
some features and observation were discussed in the final section. 
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2 The Model 


The usually input port was identified as “a-mode”, with an annihilation operator “a” (*'’ was 
omitted for simplicity). The usually unused port was called “b-mode”, with an annihilation 
operator “6” was omitted, also). Both satisfy the canonical commutation relations, say, 
[a, a*] = 1 and [6, tf] = 1. Beam splitter played the role of a mixer here. It combined both inputs 
from a-mode and b-mode then the mixture was sent into two arms as ax-mode and a 2 -mode. 
Therefore, 

a = (ai + a 2 )/\/2, b = (ai - a 2 )/\/2. (1) 

A phase was chosen. Nevertheless, it losts no generality. The reflection coefficient of arm one is 

Ri = exp (i arg r + 2 ikz\) (2) 

where “arg r” is a constant phase (real) and 2 z\ is the optical path length of the first arm. k is 
wave vector as usual. Since z is a hermitian operator, R is a unitary operator. The annihilation 
operator on one of the output ports, which was named “d-mode” as in Ref. [1] [2] , would be a linear 
combination of those from both arms. The other output port was called “c-mode” . Therefore, 

d = (Ridi + R 2 a 2 )/\/2 and c = e , ^’ c (R 1 a 1 — R 2 a 2 )/\/2 (3) 

where the phase term e l ^ c was kept for generality. Energy conservation was fulfilled. 

To treat photon shot noise and the fluctuation or radiation pressure separately was criticized 
by Marx [7] on the ground that it seems to assume some knowledge of the routes through the 
interferometer followed by individual photons, which is contrary to our understanding of quantum 
mechanics today. Loudon proposed a unified calculation [1] by introducing a coupling constant C. 
Ni pointed out that the position of the mirror itself should be a quantum-mechanical operator [2], 
It’ll contribute its intrinsic uncertainty to the total quantum uncertainty of the position of mirrors. 
In the low intensity limit, it was shown that the total uncertainty can be expressed as the sum 
of all three noise sources. Situation gets complicated at high intensity. The independence and 
correlation at different intensities among these noise moments provide ways to monitor and control 
a macroscopic quantum-mechanical object [3]. 

The Loudon-Ni model can be rephrased as the following: 

kz\ = kzu — C[a^a\ and fcz 2 = &2 2 t — (4) 

where C\ is the coupling constant which might be different from each mirror. The prime was to 
differ our notation from the previous one. There might be a factor of 2 difference, zu corresponds 
to the position operator of the first mirror without including the coupling effect, z is the mirror 
position operator we measured finally. 

The coupling constant C can be estimated as below. Suppose the mirror was hanged as a 
simple pendulum with mass M and length l. Its restoring force would be Ax-Mg/7 where Ax is 
a. small displacement and g is the gravitational acceleration. Each photon suffers a momentum 
change 2 tik after been reflected back from the mirror. On balance we got C = 2hk 2 l/(Mg). 

Photon detector usually has its own quantum efficiency, denoted by £, which was assumed 
identical for both c-mode and d-mode. The measured photon intensity would, therefore, differ 
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from (d)d) and {c*c) by a factor £. We’ve considered two detection schemes as in Ref. [8]. The 
first is direct detection: 

(m)dir = (m d ) = £{<fld). (5) 

The other is difference detection, which is the difference between the two output ports: 

(m ) diff = <m c ) - {m d ) = £((c f c) - (d'd)). (6) 

Only the variances of difference detection were presented in this article: 

(Am)d lff = £ 2 ((c f c f cc) + {dldldd) - 2{c t d t cd) - {c'c) 2 - {did) 2 + 2{c t c)(d t d)) 

+ ^({c t c) + (d t d)) (7) 

or, equivalently, in its expansion form: 

(Am)j iff = C 2 [2(a t 1 ai4a 2 ) - ((alRlR 2 a 2 ) + (4R2R-i a i)) 2 

+ (4R{a}R|R.2a2R2n2) + (a2^2 a 2^2R'i a i^-i a i)] + £(( a i a i) + ( a 2 a 2))- (8) 

These are what we want to calculate with our various inputs. It would be more complicated and 
model-dependent when considering photons with different frequencies. 

3 The Solutions 

First, a solution for coherent state - squeezed state inputs was solved. We’ll have to deal with an 
expectation value, (exp [A(oda 4- b^b) -I- B(tfa + od6)]), where the state vector | ) a is in coherent 
state |ct), | ) b is in general squeezed state | /?,(). C = se ie , where s is squeezing factor and 6 is 
squeezing angle. The coherent parameters a and 0 are complex numbers with their phases (fi a 
and (j)p respectively. 

Since coherent states |a) is a’s eigenstates, it is reasonable to reorder those operators as 

exp [ J 4(a t a + 6*6) + R(6*a + a*6)] = exp (Ua'b) exp(Fa t a) exp (T6*6) exp {Zb 1 a). (9) 

To get the coefficients U, V , Y, Z, it would be much easier to use faithful matrix representations of 
those four operators. Suppose X n , X 22 , X l2 , X 2i are their corresponding matrices, which satisfy 
the same commutators as ala, 6*6, a*6 and 6* a do. It is not difficult to find a set of faithful matrices 
(2 x 2) which have the same relations. The operators equation becomes a matrix equation after 
this substitution. Solving this matrix equation we got 

U = Z = tanh B, V = A - In (cosh B) and Y = A + ln(cosh B). (10) 

After reordering, the calculation of the expectation value on a-mode (coherent state input) 
can be carried out. (a| exp(Wa)|a) was given in Ref. [2]. What left would be a calculation 
on b-mode (squeezed state input), which looks like {fi,(\exp(Ua*b) exp(Yb ] b) exp(Uatf)\j3,(). A 
squeezed state can be expressed as a vacuum state operated by a squeezing operator S(s,6) and 
a displacement operator D(a, ct) where 

D(b, 3) = exp (06* - 0*b) and S(s, 9) = exp[{s/2)(e- 2l0 b 2 - e 2 '% n )\. (11) 
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Substitute this definition of squeezed state into b-mode, with some algebra, we may express its 
expectation value as the vacuum expectation value of a product of seperating terms 

(0X\e Ua ' b e Yb \ Uab ' |/3,C) = {0|e^( 62 / 2 >e flfrt e C6 e y ( 6t6+1 / 2 ) e D6t e Efc e F ( 6,2 / 2 )|0) x (a number). (12) 

To solve this, we turned those operators into their normal ordering. Those operators form a Lie 
algebra. With corresponding commutators and their structure constants it is possible to find a 
faithful matrix representation [6]. Therefore, the operator equation could be reordered as 

(0| exp[A(b 2 /2 )] expfBft^ exp[C7>] exp[Y (tfb 4- 1/2)] exp[£)^] exp[F6] exp[F(b t2 /2)]|0) 


= (0| exp[J(6 t2 /2) + Lb\ exp [M(tfb + 1/2)] exp[N(b 2 /2) + Pb] exp[Q]|0) 

= exp(M/2) exp(Q) (13) 

and, its corresponding matrix equation can be solved easily. We got 

g _M _ e -y(i _ AFe 2Y ) (14) 

Q = (e M /2)(2CD + 2CEF + 2ABD + 2ABEF + AB 2 e~ Y 

+ E 2 Fe- Y + C 2 Fe Y + AD 2 e Y + 2ABCFe Y + 2ADEFe Y ). (15) 

It is now straightforward to evaluate the uncertainty of photon measurements. The expectation 
value of the photon number of d-mode is 

(ftd) — (l/2)((n) + |/?| 2 + | | 2 ) + (Ho|ho|/2){cos2<//) (16) 


where |fi 0 | is roughly proportional to the input intensity and Eq is an exponential factor which 
would be discussed later. p s is squeezing related and ft is essentially the difference in optical path 
length with additional terms. 


ho — \ho\e' H ° — (n) + 


ft = (1/2) [H 0 + 2 k(z 2i -hi)- Im(Q)] 

ag’ - a*f - \p s \ 2 e Y ' fg’ 


I Ps 


\ 2 e 2Y> 


(p 2 c -\Ps\ 2 e 2Y ') 2 


(17) 

(18) 


Eo = ( Pc - \Ps\ 2 e 2Y ) * exp ( i -Vpe 2 ^ ^”)^ 1 + ~ * r ^ 1 + e ^ COS ^ “ 2 &»)) 


|r| 2 e 


+ 


m 2 (l - \r\ 2 e Y + |r|(l - e y ) cos(0 - 2*,))]} (19) 


Im(Q) = 


-iUe Y 


[(a*p + a/T)(l - |r| 2 e y ) + (a*/TT + a 0T* )(1 - e Y )\ 


1 - \Y\ 2 e 2Y ' 

f = (p 2 - \p s \ 2 e Y ')P + p c p s ( 1 - e Y ')P* + \p s \ 2 e 2Y 'U'a - p c p s e Y 'U'a* 
g‘ ' = (Pi ~ \Ps\ 2 e Y ')0* + p c p*( 1 - e Y ')p + \p s \ 2 e 2Y 'U'a* - p c p*e Y 'U'a 
U' = i tan EC', Y' = -iAC' + ln(cos EC') 

U = -i tan(EC'), Y = iAC' + ln(cosEC') 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 
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where p c = cosh s, p s = e l 6 sinh s, T = p s /p C i AC' = C' 2 — C[ and EC" = C' 2 + C[. Similarly, 

(c f c = — Eg |/io 1 (cos 20') (25) 


for difference detection. Basically, the Michelson interferometer is a transducer which turns a 
change of the arm length into a change of light intensity. The measured intensity (m) and its 
variance can be transformed back to the uncertainty of arm length, or more precisely, the difference 
of the positions of two mirrors. We may write down this uncertainty as 


( AZi a i a i )diff V , (Am& ff /(2^£'„|/.„||(sin2«|) 


(26) 


where (Am)J jf[ was given by Eq. (7) or, more explicitly, by Eq. (8). This is just the inverse of 
signal to noise ratio. The final result is 


(A Z tota i) 4- ft — 


((n) + 1/?| 2 ) 2 + 3|p s | 4 + H 2 - 4(n)|/?| 2 cos 2 (0 o - 0,) 

8A; 2 Eq | h 0 1 2 (sin 20') 2 

|/?] 2 [2|ps| 2 - pc\ps\ cos(6> - 24>p)\ + { n)pc\ps\ cos(6> - 20 Q ) 
+ 4k, 2 El |h 0 | 2 (sin 20') 2 


+ 


E" | h 2 1 (cos [40' + Hi - 2Hp + Im(Q") + 2Im(Q)]) 


(cos 20') 


8/c 2 Eq |h 0 | 2 (sin 20') 2 
A2 + (n) + \P\ 2 + \ps\‘ 


(27) 


4h 2 (sin 20') 2 4^ h 2 Eq | h 0 1 2 (sin 20') 2 


for difference detection. Where ( n ) = |a| 2 is the intensity of the input coherent state on a-mode. 
h 2 is shorthand notations of complicated modification on intensity square, H 2 is its phase. E" is 
another exponential factor which decreases the interference terms in A Ztotai- I m (Q) an d Im(Q' ) 
are imaginary parts of Q and Q", which came from solving the matrix equation Eq. (13) or its 
similar version. All of them can be evaluated exactly as follow. 


h 2 — \h 2 \e 


iH 2 


[(«) 


+ 


ag" - a*f" - 2|p s | 2 e 


2 „Y" 


f"g" 


2 , p 2 c \ps\ 2 - 2\p* 


P 2 c - I P. 


\2„2Y" 


(p 2 - Ip. 


PcPs 


n,2 — 


Ip. 


I2-2V 


7r(a ~ 


/" 


Pc - Ip. 


\2 0 2Y 


77? 


PcPs 


\2 e 2Y"y 


] + 


W r " 


(pi - \Ps\ 2 e 2Y "Y 


Pi - |p,| 2 e 2 ^' r(a * + Pi 


\Ps\ 2 e 2Y 


77? 


(28) 


U" and Y" are similar to U' and Y' but with AC' and EC" replaced by 2AC' and 2EC". /", g ", 
E" and Im(Q") are similar to /, g, E 0 and Im(Q) except Y(Y') and U{U') were replaced by Y" 


and U" . In our calculation, it was assumed that C[ = C 2 = C'. 

We now turn to the photon number states input. Photon number state is a quantum state 
without classical correspondence. Its second order coherence is minimum such that its number 
variance vanished. 

Suppose the input state is |n Q ) a 0 |nf,)fc with |n a ) a = (o^) n “/\/^J|0) |?b>)6 = (^"'’/■v/wjJlO) 

where n a and n b are real numbers. With a different reordering from Eq. (9) and the assumption 
AC' = 0, we have 


(d'd) 


^(n a + n b ) + ^(cos [2 k{z 2i - ii l )])(cos2C , ) nb n “ 

x {n a cos 2C' 2 E 1 (1 + n b , 1 — n a , 1; sin 2 2C') - -n a , 1; sin 2 2C')} (29) 
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where 2 F X (a, 6, c; 2) is hypergeometric function. In perfect detection, that is, f = 1, the variance 
of the photon detection becomes 


{(Am) 2 ) = 


^n a (n a + 1) + in b (n 6 + 1) - ](cos[2k(z 2l - £ li )]) 2 (cos2C') 2(nfc - n “ ) 

o o 4 


x[ra a cos2C" 2 Fi(l + n b , 1 — n a , 1; sin 2 2 C") orv 2 ^ 1 ( n< ” ~ n<1 ’ s ^ n2 ^C')\ 


v\i2 


cos 2 C' 


+i(cos [4fc(i 2l - 2i i )])(cos4C") nt 

O 


-n a 


X 


n a (n a — 1) cos 2 4C' 2 Fi(1 + 74, 2 — ra a , 1; sin 2 4C') 

+ ra 6 (rat - l)sec 2 4C , 2 F 1 (rat - l,-ra a ,l;sin 2 4C') 

- 4ra a rat 2 F 1 (rat, 1 - ra Q , 1; sin 2 4C") 

- ra a ra 6 (ra a - l)(ra fc - 1) sin 2 4C" 2 F!(1 + 74 , 2 - ra a , 3; sin 2 4C")] . 


(30) 


There are relations between a hypergeometric function and its contiguous functions. Further 
simplification is possible. 


4 Discussion 

A low intensity limit can be obtained easily [1] while b-mode was in vacuum state. In short, 
(AZtotai ) 2 > C'/k 2 . We have searched in a limited parameter space and found no violation of 
this inequality (within error). Though, it is not a proof of that limit. At very high intensity, the 
reflectance R has little effect. The noise behavior can be explained by the mixing of two input 
states. High order moments are needed to characterize such a superposition. On the other hand, 
that noise can be eliminated by setting two input states at nearly the same intensity. Nevertheless, 
there is still an exponential factor Fo in signal (cf. Eq. (25) and Eq. (26)). It doesn’t show up in 
classical solution and it always decreases the signal. We left further discussion to another work. 
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